Thermal Effects on the Stability of Excited Atoms in Cavities 
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An atom, coupled linearly to an environment, is considered in a harmonic approximation in 
thermal equilibrium inside a cavity. The environment is modeled by an infinite set of harmonic 
oscillators. We employ the notion of dressed states to investigate the time evolution of the atom 
initially in the first excited level. In a very large cavity (free space) for a long elapsed time, the 
atom decays and the value of its occupation number is the physically expected one at a given 
temperature. For a small cavity the excited atom never completely decays and the stability rate 
depends on temperature. 



I. INTRODUCTION 

Inhibition of spontaneous emission by confined atoms 
is a well-known phenomenon, currently related to the 
dipole orientation with respect to parallel mirrors or to 
the relation between the size of the confining, device and 
the emission wavelength (see for instance [H,0| and other 
references therein). The theoretical understanding of 
these and other effects in atomic physics on perturbativc 
grounds requires the calculation of very high-order terms 
in perturbative series, that makes the standard Feyn- 
man diagram technique practically unreliable. This has 
lead to trials of treating non-perturbativelly such kind 
of systems using the semi-quantitative idea of a dressed 
atom [3|. However serious difficulties, due to nonlinearity, 
are present to get rigorous results in these approaches. 
A way to circumvent these mathematical difficulties is to 
assume that under certain conditions the coupled atom- 
electromagnetic field system may be approximated by a 
system composed of a harmonic oscillator coupled lin- 
early to the field modes through some effective coupling 
constant g. This is the case for linear response theory in 
QED, where the electric dipole interaction gives the lead- 
ing contribution to the radiation process [J, [5| . Although 
a linear model is a simple theory, it permits a better un- 
derstanding of the need for non-perturbative analytical 
treatment of coupled systems. This is the basic prob- 
lem underlying the idea of dressed quantum mechanical 
operators. 
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The perturbative treatment of interacting systems is 
carried out by considering bare, non-interacting fields. 
The interaction is taken into account order by order in 
powers of the coupling constant. However there are sit- 
uations where the use of perturbation theory is not reli- 
able, as in the low energy domain of quantum chromody- 
namics and resonant effects in atomic physics, associated 
with the coupling of atoms with strong radiofrequency 
fields 0. 

The idea of a bare particle associated to a bare matter 
field is actually an artifact of perturbation theory. A 
charged physical particle is always coupled to the gauge 
field, i.e, it is always "dressed" by a cloud of quanta of 
the gauge field (photons, in the case of electrodynamics). 
In a simplified model for a radiating atom, a way to treat 
directly dressed objects, has been introduced. This is the 
method of dressed states and dressed coordinates 0| that 
has been employed in several cases [7MTTI . 

In this paper we generalize the zero-temperature for- 
mulation, dealing non-perturbatively with the inhibition 
of spontaneous emission [7|,|8j, to finite temperature. The 
objective is to consider the stability of confined excited 
atoms in an environment at finite temperature. 



II. THE MODEL 

We start by considering a bare atom approximated 
by a harmonic oscillator described by the bare coordi- 
nate and momentum qo,po respectively, having bare fre- 
quency wo , linearly coupled to a set of N other harmonic 
oscillators (the environment) described by bare coordi- 
nate and momenta qk,Pk respectively, with frequencies 
u>k, k = 1,2, . . . , N. The limit N — > 00 will be taken 
later. A model of this type, describing a linear coupling 
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of a particle with an environment, has been used for years 
in several situations, for instance to study the quantum 
Brownian motion of a particle with the path-integral for- 
malism [l2T - fl5l ]. The whole system is supposed to reside 
inside a spherical cavity of radius R in thermal equilib- 
rium with the environment, at a temperature T = /3 _1 
(fcg, the Boltzmann constant is taken equal to 1). 

The Hamiltonian for such a system is written in the 
form, 
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where the c^-'s are coupling constants. In the limit 
N — » oo, we recover the case of an atom coupled to the 
environment, after redefining divergent quantities, in a 
manner analogous to mass renormalization in field theo- 
ries. 

The Hamiltonian ([T]) is transformed to the principal 
axis by means of a point transformation, 
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where \i — (0, {fc}) k = 1,2, ...,N, performed by an or- 
thonormal matrix T — (tV). The subscripts /i — and 
fi = k refer respectively to the atom and the harmonic 
modes of the reservoir and r refers to the normal modes. 
In terms of normal momenta and coordinates, the trans- 
formed Hamiltonian reads 
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where the f2 r 's are the normal frequencies corresponding 
to the stable collective oscillation modes of the coupled 
system. It can be shown [6] that 
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The sum in Eq. (HJ diverges for N — > oo. This makes 
the equation meaningless, unless a renormalization pro- 
cedure, analogous to mass renormalization in field theo- 
ries, is implemented 16]. Adding and subtracting a term 
r] 2 Vl 2 in the numerators of the right hand side in Eq. (|4]) 
we have 



j 2 - n 2 r = r] 2 n 2 r 
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where we define the renormalized frequency 
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We find that the addition of a counterterm — Suj q$ in the 
Hamiltonian Eq. (TTJ compensates the divergence of Wg in 
such a way as to leave a finite, physically meaningful 
renormalized frequency uj. 
Using the formula, 
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Eq. ((7} can be rewritten as (dropping the label for the 
eigenfrcquencies) , 
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This gives an infinity of solutions. The spectrum of 
the collective normal modes is denoted by tt r ; r = 
0, 1,2, • • •. The transformation matrix elements are @, 

rjQ r 
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Unless explicitly stated, the limit N 
in the following. 
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To correctly describe the coupling of the atom with the 
field, we take 



c fe = r)u k ; rj 
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where g is a constant with dimension of frequency. It 
measures the strength of the coupling; Auj = ixcjR is the 
interval between two neighboring frequencies of the reser- 
voir and frequencies of the field modes are given by , 
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III. DRESSED STATES 

Let us consider the eigenstates of our system, 
|riO) n>i, n.2...}, represented by the normalized eigenf unc- 
tions in terms of the normal coordinates {Q r }, 
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where H ns stands for the n s -th Hermite polynomial and 
Tq is the normalized ground state eigenfunction, 
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We introduce dressed coordinates q' and {q' k } for, re- 
spectively, the dressed atom and the dressed field, de- 
fined by, 



(14) 



where u)^ = {ui, uj k }- In terms of dressed coordinates, we 
define for the time r = 0, the dressed states, \kq, ki, k.2---) 
by means of the complete orthonormal functions 
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where q' = {q' Q , q' k }, d> M = {ui, uik}- Notice that the 
ground state, To, in the above equation is the same as 
in Eq. (| 12[) . The invariance of the ground state is due to 
our definition of dressed coordinates given by Eq. (1141) . 
In fact, we get the normal coordinates Q r in terms of the 
dressed ones from Eq. (fT4")) . Replacing them in Eq. (TIB"]) 
we find that the ground state in terms of the dressed 
coordinates has the form 
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Each function ip K , oK , 1 ...{q') describes a state in which the 
dressed oscillator q'^ is in its n^-th excited state. 

It is worthwhile to note that our dressed coordinates 
are new objects, different from both the bare coordi- 
nates, q, and the normal coordinates Q. In particular, the 
dressed states, although being collective objects, should 
not be confused with the eigenstates given by Eq. (TT2"1) . 
While the eigenstates <f> are stable, all the dressed states ip 
are unstable, except for the ground state Tq . The impor- 
tant idea is that the dressed states are physically mean- 
ingful. 

In this framework, we write the physical states in terms 
of dressed annihilation and creation operators a'^ and a'} 
defined in terms of dressed coordinates and momenta in 
the usual way, 
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Then the initial dressed density operator corresponding 
to the thermal bath is given by 
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with Z'p = Y\ k Zp' being the partition function of the 



dressed reservoir, where 
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The system evolves with time (r). The time-dependent 
dressed occupation numbers are defined as 

(the prime is to clearly distinguish the dressed quanti- 
ties from the bare ones), where p' is the density opera- 
tor for the dressed atom; a'(r) and a^(r) are the time- 
dependent creation and annihilation operators. 

The time evolution of the dressed annihilation operator 
is given by, 
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and a similar equation for aj|(r). We solve this equation 
with the initial condition at time r = 0, 




which, in terms of bare coordinates, becomes 
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We assume a solution for a'(r) of the type 
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Using Eq. (fTJ) we find, 
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The initial conditions for B' (r) and B' (r) are obtained 
by setting t = in Eq. (|24[) and comparing with Eq. 
then 



r=0 



j-r u-r 



KM 



Km = EVv^ 



r=0 



(26) 



(27) 



Using these initial conditions and the orthonormality of 
the matrix {t 1 ^} we obtain a'£ — 0, b'£ — itTj s/Wjj.. Re- 
placing these values for a? and in Eq. (j25]l we get 
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where 



(30) 



with fi , v = 0, {fc}, k = 1, 2, • • •. 

This leads to the time evolution equation for the 
dressed occupation number of the atom [T{J [2(| , 

oo 

n , Q (T) = \foo(r)\ 2 n' + Y / \fok(r)\ 2 n / k , (31) 

fe=i 

where n' stands for the occupation number at r = 0. 

IV. THERMAL EFFECTS IN A SMALL CAVITY 

In this section we consider the weak coupling regime, 
defined by 

g = u>a,, (32) 

where a is the fine-structure constant. 

With 77 = ^2g-Kc/R and defining the dimensionless 
parameter 
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Eq. (fTU)) becomes, 
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Let us consider the right hand side of Eq. (|34l) such 
that 
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In the weak coupling regime, this corresponds to a value 
of 5, 5 > 5.3 x 10~ 5 . For a frequency w = 4.0 x 10 14 /s 
(in the visible red) this gives a condition on the cavity 
size of R > 1.7 x 10~ 8 to. Then the general behavior of 
the solutions of Eq. (l34|) is illustrated in Fig(p}. We find 
that all but one of the eigenfrequencies are very close to 
the frequencies of the field modes, lj^, given by Eq. ([B]). 
Then we label solutions for the eigenfrequencies £l r as, 
il , ilk, k = 1,2, • • •. The solutions Slfc of Eq. (J34J) are, 
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Since every et is much smaller than 1, Eq. (|37[) can be 
linearized in e/c, giving, 
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FIG. 1: Solutions of Eq. (|34[) . for cavities satisfying the con- 
dition Eq. (|35[) ; the asymptotes of the cotangent curve corre- 
spond to the frequencies of the field modes ujk- 



The eigenfrequencies, f2&, are obtained by solving 
Eqs. (|36t and (|3T |) or (|38 ]l . 

The lowest eigenfrequency, £7 , is obtained by assuming 
that it satisfies the condition tt R/c <§; 1. Inserting this 
condition in Eq. (|34[) and keeping up to quadratic terms 
in f2 we obtain the solution for the lowest eigenfrequency, 
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Consistency between Eq. (|34|) and the condition 
QqR/c <C 1 gives a condition on R, i.e R <C (c/g)X, 
with A = (tt/2) (g/w) 2 . 

Let us first determine the temperature independent 
term |/oo(t) | 2 "-o in Eq. (13T1) . considering that the dressed 
atom is initially (at r = 0) in the first excited level, that 
is n' = 1. We evaluate (i°) 2 and (i§) 2 from Eqs. (fTT]) . 
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and then using the normalization condition X)^Lo(^o) 2 



1 and C(2) = J2k=i k = 7r 2 /6, we have 
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From Eq. (|30[) . using the de Moivre formula, e = cos 6 + 
i sin 9, we have 

00 

r,s— 

Let us assume that the thermal bath is at zero tempera- 
ture, i.e. all the modes of the reservoir are in the ground 
state, n' k — for all values of k. Taking the above ap- 
proximations for t§ and ig, we get from Eq. (|3"Tj) . the 
zero-temperature time evolution of the occupation num- 
ber of the dressed atom initially in the first excited level, 
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This is an oscillating function which has a minimum value 
Min(|/ 00 (*)| 2 ). Taking both cosine functions in Eq. fl33J) 
equal to — 1, we get a lower bound for Min(|/ 00 (i)| 2 ) 
given, up to first order in S, by 
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As an example we consider that the atom in the first 
excited state has an emission frequency ui ~ 4 x 10 /s, 
in the visible red, and we take the radius of the confining 
cavity R 10~ 6 m. With these data we get F(S) « 0.99, 
that is a probability of 99% at zero temperature, that it 
will almost never decay. This shows the high stability of 
the system, which is confirmed by experiment (T3, EH ■ 

In order to take into account the temperature effects, 
we must consider the second term in Eq. (1311) , that is we 
must evaluate the quantity 
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This is carried out by using the matrix elements obtained 
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from Eq. ([lip and the formulas for eigenfrequencies in a 
small cavity. 

It is assumed that the thermal distribution of the oc- 
cupation numbers of the field modes in the cavity follow 
the Bose-Einstcin distribution, 
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This can be justified in the following way: in the case of 
an arbitrarily large cavity, the dressed field modes coin- 
cide with the bare ones |6j , which in the limit of vanishing 
coupling makes this approximation exact. Strictly speak- 
ing this is not the case for a finite cavity. Nevertheless, 
in many situations this approximation is acceptable in 
the weak coupling regime. For instance a cavity of ra- 
dius R « 10 _6 m is ~ 10 4 times larger than the size of 
a hydrogen atom (the Bohr radius). In such a case the 
atom "sees" the cavity to be a very large one and the 
approximation is justified. Then from Eqs. (H51 and (|3ip . 
we get the time evolution of the temperature dependent 
occupation number for the atom, 
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where |/oo(t)| 2 is given by Eq. (|43|) . The matrix ele- 
ments iS and t l k in the above formulas are evaluated from 
Eqs. (HO), (ES) and O to be, 
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k 2 g 2 - tt 2 6 2 ' k k 2 -{l + ei) 2 I ' 
The occupation number n/ (r, j3) is an oscillating function 



which has a minimum value, Min[no(r, that depends 
on the temperature We can obtain a lower bound, 
n' (/3), for this minimum, such that Mm[n' Q (T, (3)} > 
n' (f3), by taking both cosine functions in Eq. (|4"7|) equal 
to -1 
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Numerical calculation of Eqs. (|47[) and (|49p will describe how the time evolution of the occupation number and 
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be attained in laboratory for excited atoms. In fact it 
is lower than the ionization temperature of 13.6 eV for 
the hydrogen atom, and still much lower than the nu- 
clear fusion temperature of « 10 8 K (~ 8AKeV). We 
find that the average occupation number at this temper- 
ature (T = 10 5 K ) is about 3.5 times higher than the zero 
temperature value n' Q (T,T = 0) w 1. At room tempera- 
ture the occupation number will remain very close to the 
zero-temperature value, as is shown in Fig. ([2]). There- 
fore we find that as the temperature is raised, both the 
amplitude of oscillation of the occupation number and 
its minimum, grow with respect to the zero-temperature 
values. 



FIG. 3: Time evolution of the thermal dependent occupation 
number n' (T,f3) for a temperature T — 10 K; time is in 
seconds. 



the stability of the excited atom are affected by heating. 
We take for the plots n' = 1, that is the atom initially 
in the first excited level. 

In Fig. @ and in Fig. © the time evolution of the 
temperature dependent occupation number 71q(t, fi) is 
plotted for some representative values of the emission 
frequency and the temperature. In Fig. (j4]) the lower 
bound for its minimum, n' (/3), is plotted as a function 
of temperature. We find from these figures that raising 
temperature increases the amplitude of oscillation of the 
occupation number and that its minimum lower bound, 
n' (/3), also grows with temperature. For a given emission 
frequency ui = 4.0 x 10 14 /s, the increase of the amplitude 
of oscillation of ti' q (t,(3) and of its lower bound n' (j3) 
are negligible for room temperatures, they are significant 
for high laboratory temperatures. In Fig. ([3]) n' Q (T,f3) is 
plotted for T = 10 5 K (~ 8.4 eV); although this temper- 
ature is very high for everyday life, it can, in principle, 



V. CONCLUDING REMARKS 

At zero temperature the dressed atom, initially in the 
first or higher excited state, can only decay, since all field 
modes are in the ground state. It is inhibited from de- 
caying by confinement in a cavity of small size. However 
at finite temperature, the field modes in the cavity can 
be in excited states with a finite probability given by the 
Bose-Einstein distribution function. As a consequence 
the dressed atom can exchange quanta from the field. 
This means that the thermal occupation number of ex- 
cited states of the atom increases with temperature. In 
other words, as an effect of heating the atom will be in 
a higher excited state which may be able to decay. How- 
ever the decay is inhibited by the confining geometry. 
The results presented above give sufficient proof of these 
ideas. 

This behavior is also to be contrasted with the situa- 
tion of an arbitrarily large cavity (free space) described 
in [20| . In that case, for long times the dressed oc- 
cupation number of the atom approaches smoothly to an 
asymptotic value which is nearly the one obtained from 
the Bose distribution at the equilibrium temperature of 
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the reservoir. Taking the same value as before for u>, this 
value is n' (t — > oo, f3; ui, R — > oo) = l/(e 14 - 1) » 0. In 
that case the growth of the Bose-Einstein weight factor 
due to raising temperature is compensated by the low- 
ering due to larger and larger values of R, leading to an 
equilibrium occupation number. 
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